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Investigation of the Fusibility of 
Twp>r«turt Prof Hint Optical Diagnostics 
in the SSME Full Fre-Burner 




Results of sn snslyticsl invest igat ion performed by United Technologies 
Research Center (UTRC) under Contract HAS8-34744, sponsored by the George C. 
Marshall Space Flight Center to determine the feasibility of temperature 
profiling in the space shuttle main engine (SSME) fuel pre-burner are presented. 
In this application it is desirable to measure temperature in the pre-burner 
combustor with a remote, non- intrusive optical technique. Several techniques 
using laser excitation have been examined with a consideration of the constraints 
imposed by optical access in the fuel pre-burner and the problems associated with 
operation near the functioning space shuttle engine. The potential performance 
of practical diagnostic systems based on spontaneous Raman backscattering, laser- 
induced fluorescence, and coherent anti-Stokes Raman spectroscopy have been 
compared analytically. A system using collection of spontaneous Raman back- 
scattering excited by a remotely- located 5-10 Hatt laser propagated to the SSME 
through a small diameter optical fiber was selected as the best approach. 
Difficulties normally associated with Raman scattering: weak signal strength and 
interference due to background radiation are not expected to be problematic due 
to the very high density in this application, and the low fleam luminosity 
expected in the fuel-rich hydrogen-oxygen flame. The components of a Raman 
system have been examined critically and the limitations, if any, that they 
impose have been considered. Theoretical estimates of the performance of the 
Raman system indicate that it should be capable of achieving system design goals 
with respect to spatial, temporal, and temperature resolution with minor 
tradeoffs in overall performance. 
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INTRODUCTION 


The Space Shuttle Main Engine (SSME) ia one of the most critical components 
of probably the most ambitious engineering endeavor of the past decade. At full 
power it develops 512,000 lbs of thrust, and can be throttled frosi 65 percent of 
the 470,000 lbs of rated thrust to 109 percent. The SSME is designed for a life 
of 7.5 hrs with 55 starts, with some servicing of line replaceable units. Major 
technical hurdles have had to be overcome to develop this engine. Early in the 
development of the engine, problems were experienced with the failure of bearings 
in the oxidizer turbopunp. Turbine blade life has been an area of concern and is 
the subject of an on-going upgrade of engine performance in terms of extending 
the usable life. This report indirectly addresses this problem. Turbine blade 
life is thought to be limited by fatigue possibly caused by nonuniform tempera- 
ture profiles in the pre-burner. The development of modifications leading to 
more uniform temperature profiles would be greatly facilitated by a diagnostic 
capable of measuring temperatures in the engine. The objective of the investiga- 
tion reported herein is to determine the best optical techniques for temperature 
profiling in the SSME pre-burner, and to determine on theoretical grounds the 
probability of success. 

Optical diagnostics have been under development since the advent of high 
power and tunable lasers. A number of techniques have evolved to make measure- 
ments in a wide variety of combustion systems, including such practical environ- 
ments as combustion research tunnels, internal combustion engines, large 
industrial furnaces and after-burning jet engine exhausts. Special, difficult 
problems are encountered in the environment of the SSME. These include: engine 

motion; fog, mist and ice formation near and on the engine; and severe vibration 
and acoustic noise levels in the vicinity of the engine. Optical measurements 
are also hampered by existing optical access into the pre-burner. Several 
optical techniques have been examined in light of these restrictions. These are 
laser-induced fluorescence, spontaneous Raman, and coherent nonlinear techniques, 
such as coherent anti-Stokes Raman spectroscopy (CARS). The potential perfor- 
mance of these optical techniques has been predicted and compared, based on 
properties of the flow in the pre-burner combustor and theoretical predictions 
and experience in optical combustion diagnostics. 

Spontaneous Rasan scattering has been selected as the best available 
technique based on the restrictions posed by the geometry and environment of the 
SS>E pre-burner and the predicted favorable performance of Raman diagnostics. 
Spontaneous Raman requires only single port access, is capable of adequate 
performance in terms of program goals, and is amenable to implementation in an 
optical fiber configuration. This latter point is very important because the 
severe environment of the SSME engine dictates taking steps to maintain measure- 
ment integrity. An optical fiber should be insensitive to engine motion, 
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condensation on and around tha angina* and tha noise and vibration levels near 
the angina. Coherent optical techniques would be capable of aeasureaents only 
along a line joining tha two existing access ports in tha coabustor. Fluores- 
cence aeasureaents are not favorable due to a lack of adequate concentrations of 
species that are aaenable to laser excitation. The introduction of a seed 
species is probleaatic and the existence of other usable techniques argues 
against this node of operation. 

The theory of Raaan aeasureaents in hydrogen* the aajor constituent in the 
S8MS pre-burner poet coabustion flow* has been exaained and the aystea eleaents 
coaprising the diagnostics have been identified. In particular the use of 
optical fibers in this application has been closely studied. The levels of 
performance predicted considering practical system limitations compare favorably 
with the design performance criteria established by NASA. Probleas which have 
not been studied are those which cannot be effectively studied on theoretical 
grounds alone. In particular* one potential problea which asy hamper optical 
diagnostics is associated with the transaission of laser beams through the SSME 
pre-burner due to possible severe coabustion turbulence or two-phase flow 
problems. In this regard the selected technique appears to be the most 
reasonable approach. 

The results of the feasibility study are described in the following sections 
of this report. The next section describes the relevant optical techniques for 
SSME pre-burner optical diagnostics. The section opens with a description of the 
fuel pre-burner especially as regards optical probing. Then the diagnostic 
techniques are described and the performance of each is coapared. The next 
section focuses on the use of Raman scattering techniques. The section starts 
with a review of the principles of Raaan temperature aeasureaents. The component 
parts of the Raman system are examined next. Particular attention is psid to the 
use of optical fibers in this section. Finally in this section* the performance 
of an optical diagnostic system based on Raman diagnostics is exaained in the 
light of limitations of system elements just studied. Conclusions drawn from the 
study sre then presented. References are contained in a separate section near 
the end of the report. Appendix A contains a compilation of Raman optical 
intensity calculations. 
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OPTICAL DIAGNOSTIC TECHNIQUES FOR TEMPERATURE MEASUREMENTS 


This section describes the verious optical diagnostic techniques available 
for teaperature measurements in high pressure combustion. The 8SME fuel pre- 
burner is described first, since its characteristics influence the choice of 
optical techniques. After this the optical diagnostic techniques most relevant 
to this application are introduced. The perforaence of these techniques is 
enuaerated in the subsection following end the selection of best optical 
technique is justified. 

The perforaence capability end the liaitations of the selected technique, 
Raaen backseat ter ing is described in the section following this one. 


Constraints on SSME Pre-Burner Optical Meesureaents 

The fuel pre-burner reacts a large portion of the hydrogen flow with a snail 
portion of the oxygen flow to provide a hot gas streaa to drive the turbines that 
drive the fuel pump. The hot gas stream leaving the turbine is ducted to the 
aain combustion chamber where it reacts with the remainder of the oxygen flow. 

The pre-burner operates very fuel-rich so the gas temperature in the turbopump is 
relatively low, approximately 970 R (1750 R) at the nominal power level. 

Equilibrium calculations at the fuel pre-burner flow conditions at the 
nominal power level (NPL) indicate that molecular hydrogen is the dominant 
constituent in the combusted flow. Water vapor makes up the balance of the 
mixture at about 10 percent of the total flow. The atomic hydrogen, hydoxyl 
radical (OH), and hydrogen peroxide mole fractions are estimated to be very low. 
The calculated gas composition is shown in Table 1. 


TABLE 1 

FUEL PRE-BURNER CALCULATED COMPOSITION 


NPL, p ■ 5500 psia 


Species 


Mole Fraction 


H 2 

0.89407 

h 2 o 

0.10593 

H 

1.559 x 

OH 

4.709 x 

) 2 1 0 g fltid 0 j 

< io- 20 
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A cross section view of the fuel pre* burner is shown in Fig* 1« Fuel end 
oxiditer ere introduced into the pre-burner chewber by injection (not shown) 
through the fees piets. They reset end exheuet to the turbine. Within the 
chamber ere beffles end e liner, both of which ere supplied with hydrogen flow 
for cooling. Trensducer ports in the chamber wells permit opticel access to the 
combustion medium. A epeciel hole with e reinforcing bushing goes through the 
liner at the existing ports. There ere only two existing trensducer ports but 
two other ports could be added to bosses existing on the chamber. The ports ere 
located at 60, 165, 210 end 300 degrees with respect to the fuel inlet. 

Therefore, there is no optical access directly across the pre-burner, that is 
line-of-sight through either the existing, or the additional ports. The baffle 
plates also restrict opticel access as the axis of the trensducer ports inter- 
sects the lower edge of the beffles. It is desirable, however, to restrict 
diagnostic systems to the use of these ports to avoid unnecessary penetrations 
through the chamber. These considerations favor single port techniques. 

There ere other problems associated with the challenging application of 
optical diagnostics in the environment of a large rocket engine. First, the SSME 
moves in the test stand, because the engine gimbals for thrust vector control and 
because the thrust load deflects the test stand slightly. These motions would 
preclude free transmission of a beam to the engine due the critical alignment 
usually required of optical components. Another factor that makes transmission 
of the laser beam to the engine difficult is the fact that prior to and during 
firing, the engine is surrounded by a vapor mist due to the condensation of 
atmospheric water by the cold cryogenic propellant lines. This mist would 
scatter a laser beam severely, thereby requiring the beam to be transported in a 
special tube purged with dry nitrogen, or an optical fiber. During firing and 
afterwards, frost forms on the engine and would be expected to form on attached 
components unless measures are taken to prevent it, such as by heating with warm 
nitrogen. Finally, very high vibration and noise levels are experienced at the 
engine, and somewhat lower levels prevail in instrumentation bays located in 
nearby hardened core of the engine test stand. 

The optical properties of the flow within the pre-burner is unknown from an 
empirical point of view. From known flowrates it can be estimated that the Mach 
number is 0.005, so that compressibility effects are expected to be small. 
Monuniformitics in the refractive index resulting from incomplete combustion and 
possible two-phase flow may be problematic. The point of measurement is close to 
the point where coolant hydrogen is injected from the baffles. The magnitude of 
these possible problems can only be assessed in hot gas tests. 


SSME PRE-BURNER CROSS SECTION 
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Optical Techniques for Teaperature Measurements 

The proliferation of high power and tunable lasers has encouraged the 
development of a number of optical methods for use in hostile flasM environments 
(Refs. 1-3). Chief among the reasons for these developments is the desirability 
of non- perturbing measurement methods. The various optical techniques 
essentially probe individual excitation levels within the chemical species. The 
strengths of response of the various molecular levels can be related through 
analysis to populations in the flame from which the temperature can be deduced. 

A number of optical techniques have been developed. Among these spontaneous 
Raman, fluorescence and nonlinear optical techniques such as CARS and SRG have 
received the greatest attention recently for combustion applications. The 
advantages and disadvantages of these techniques for SSME pre-burner optical 
diagnostics are discussed in the following. 

Fluorescence 


Fluorescence (Refs. 1, 3, 4) is the emission of light from an atom or 
molecule promoted to an excited state. In laser fluorescence, the excitation 
takes place by the absorption of laser radiation tuned to coincide with a 
resonance in the probed species. The temperature is determined from an analysis 
of the distribution of spontaneous emission from the excited species. The 
hydroxyl radical is the only species existing in the pre-burner post flame gases 
at equilibrium (see Table I), which is amenable to probing by existing laser 
sources. OH has absorption bands between 2400 and 4000 Angstroms. Hydrogen and 
water vapor have optical absorption bands in the vacuum ultraviolet (Ref. 1), a 
region where few lasers sre available and no windows can be used. However, the 
concentration of hydroxyl radical is very low due to the low combustion tempera- 
ture and the fuel rich conditions in the pre-burner, making it a relatively poor 
choice as a probe species. Laser fluorescence is used to measure trace species; 
therefore, the use of this radical was investigated, and is reported in the next 
subsection. 

The penalty imposed on fluorescence may be lifted by seeding the flow with a 
chemical substance which has good fluorescence properties. Seeding may be 
accomplished by mixing the compound with one of the propellants or by injecting 
it directly into the pre-burner. Both of these approaches presents problems. To 
mix the fluorescent species with the propellant requires that it be chemically 
inert with respect to the propellant and also that it be miscible with the 
cryogenic liquid phase. Injection would require a large number of injection 
ports in the injector face, which is not practical either. 
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Raman Scattering 

Kntn scattering (Refs. 5-7) is the process by which s molecule interacts 
inelastically with sn incident photon, i.e., energy exchange between the 
molecules end photons occurs. If the molecule is initially in an unexcited 
state, the interaction leaves the aolecule in an excited state and the scattered 
photon energy is lower than the incident photon energy by an amount corresponding 
to the molecular energy. Typically a nuaber of internal energy levels can be 
exr.ited, and they each scatter light at a characteristic frequency shift so that 
the temperature can be determined from the distribution of frequencies in the 
scattered radiation. Normally only a small nuaber of molecules is excited in the 
scattering process so that no perturbation of the molecular populations occurs. 
The usual problem with Raman scattering for combustion diagnostics is that the 
scattering cross section is very small and the scattering is waak, so that back- 
ground luminosity or laser-induced interferences may overwhelm the signal. Raman 
scattering also takes place into 4v steradians, so that usually only a fraction 
of the scattered light can be collected. Fortunately the SSME pre-burner is 
likely to have very low flame luminosity because a hydrogen flame emits very 
little radiation in the visible part of the spectrum. Furthermore, the high 
chamber pressure helps considerably to overcome the small Raman scattering cross 
section. 

Nonlinear Optical Techniques 

The emergence of high power lasers has promoted the development of several 
nonlinear optical techniques for combustion diagnostics. Chief among these ere 
coherent anti-Stokes Raman spectroscopy (CARS) (Refs. 8-10) and stimulated Raman 
gain (SRG) (Ref. 11). The advantage of CARS and indeed all nonlinear techniques, 
which are coherent, is that signal radiation is generated or contained in a 
highly collimated beam which can be easily collected. Furthermore, high power 
lasers usually are employed, with which the conversion of incident radiation into 
CARS can be quite efficient. The disadvantages of these techniques is the 
limited optical access existing in the pre-burner, which restricts temperature 
profiling to a chord joining the existing access ports in the pre-burner. CARS 
also requires a knowledge of linewidths including the effect of temperature, and 
collision partner on the width. 

The performance of a CARS measurement system for this application is 
enumerated in the following section for comparison to the other techniques. 
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The advantages of tha Raman, fluoraacanca and CARS ara aummariaad in 
Tab la 11. Tha Raman tachiqua looks to bn tha moat favorabla approach for SSME 
pra-burnar diagnoatica. Tha laat row of thia chart praaanta a compariaon which 
haa not baan mada ao far. Thia ia, tha aaaa of implamantation. Tha uaa of 
optical diagnoatica in tha vicinity of a rockat angina praaanta aavnral 
challangaa already daacribad. Ona advantage of Raman acattering ia that it 
appaara quite faaaibla to uaa an optical fiber to connect to an optical head 
located on one aide of tha pra-burnar. Thia would ovcrcoma many of tha 
difficultiaa aaaociatad with operation of tha angina, auch aa tha motion 
aaaociatad with thruat load and gimballing, fronting and vapor condenaation that 
occurs around the engine due to the cryogenic propellants and vibration and 
acoustic radiation. 

Spontaneous Raman scattering appears to be the clear choice for SSME pre- 
burner optical diagnoatica if signal magnitudes are adequate, and these are 
considered in the following sections. 


Compariaon of Optical Diagnostic Techniques 

This section compares numerically the potential performance of the Raman 
scattering, CARS and laser fluorescence techniques for SSME pre-burner optical 
diagnostics. The diagnostic system performance is determined by pre-burner gas 
flow conditions and geometry. The pre-burner is assumed to be operating at the 
nominal power level at which the pressure is 38 MPa (5500 pai) and the tempera- 
ture is 1000 K (1800 R). The gas composition is taken as indicated in Table 1. 
The existing access ports are 1,1 cm (0.55 inches) in diameter and the combustor 
is about 25 cm (10 inches) in diameter. The available cw optical power in the 
pre-burner is fixed at 5 Watts, which is the approximate level allowable in an 
optical fiber of diameter appropriate for these measurements. The measurement 
spatial and temporal resolution are taken as 1 cm and 10 msec, respectively. 

Raman Scattering 


The Raman scattering signal power excited by an incident laser of power 
is given by (Ref. 1) 


P r P 

A L 





( 1 ) 


a 
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TABLE 11 

COMPARISON OF THE ADVANTAGES AND DISADVANTAGES OF 
APPLICABLE TECHNIQUES FOR 8SME OPTICAL DIAGNOSTICS 


Aspect Reman 


CARS Fluorescence 


Required geometry 

Any geometry 
satisfying 
polarisation 
requirements 

Line-of-sight 

Any geometry 

Source requirement 

Single laser 

Two lasers 

Single tunable 
laser or perhaps 
two lasers 

Collection 

Scattered into 
4s steradian 

Efficient 

collection 

possible 

Scattered into 
4s steradian 

Anslysis of date 

No theoretical 
hurdles 

Linevidth 
data required 

Quenching rates 
required 

Signal level 

Adequate 

Limited by 
approach 

Lov OH concentrations 
Other molecules 
inaccessible 
Possible seed 
required 

Implementation 

Optical fiber 
approach appears 
feasible 

Optical fiber 
not feasible 

Poor transmission of 
required UV light if 
not seeded 
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«6«ri nj is the population of the J-th rotational level, do/dQ la the 
diffarantial Reman cross sac t ion, Q is tha coll act ion solid ant la, c ( is tha 
col lac t ion afficiancy and I is tha sampling extent or aoatial resolution. The 
Renan cross section for B 2 is 1.4 x 10“*® cm 2 staradian for a incident laser 
wavelength at 4180 Angstrom (Kef. 12). Tha solid angle for a 1 cs aperture at 
13 ca, half the coabustor disaster, is 0.004 sterediana. Taking the collection 
efficiency as 10 percent and the s sapling extent as 1 ca, the calculated Raaan 
signal is 5.5 x 10* photons. This level of perforaance is shown in the next 
section to be adequate to aaet program design goals. That section further 
considers the perforaance for practical systea eleaents. 

Coherent Anti-Stokes taaan Spectroscopy 

CABS is a wave airing technique requiring two frequencies. One laser 
supplies one frequency, called the puap, while another laser is used to supply 
tha second frequency which is Stokes shifted with respect to the first. Two 
approaches are used to generate the CARS spectrua. In the first, the Stokes 
laser is scanned across the appropriate frequency interval, and the entire 
spectrua is generated successively. This approach is liaited by the tiae needed 
to scan the laser. The second approach uses a wide bandwidth Stokes dye laser to 
generate the entire aolecular spectrua at once. 

The two approaches divide into continuous-wave and pulsed aodes of operation 
because of the characteristics of dye lasers. Pulsed lasers afford gain over a 
broader spectral profile, and can have a bandwidth as broad as a few hundred 
wavenuabers, whereas cw dye lasers have bandwidths of 1 as - * or saaller. Both 
approaches have been considered for SSlfi pre-burner optical diagnostics. 

The cw dye laser approach would use two or aore dye lasers tuned to 
resonance with hydrogen aolecular transitions. The CARS signal can be estimated 
from the expression for CARS power P 3 ic terms of the laser powers (Ref. 1). 



where Sj is the puap frequency Sj is the anti-Stokes frequency, P j and P 2 are the 
puap and Stokes laser powers respectively, x i* the third order nonlinear 
susceptibility. Assuming the same puap power as in the other cases, and a dye 
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later power of 0.2 W, which it typical of available cm dye latare and a 
collection efficiency of 10 percent, the calculated CABS signal it 2 a 10 4 
phot one. It thould be noted that thit approach it cottly due the expente of 
using two dye late e. A typical commercial dye later would have a bandwidth of 
1 ct“ l and would r juire frequency ttabilitation to remain on the peak of the 
Sanaa thifted retonance, which it on the order of 0.4 ct' 1 . 

The aecond approach to CABS diagnottict utet a pul ted later to provide the 
pump frequency and to pump a broadband dye later. The utual output of a pulsed 
later it too high to intert into an optical fiber, and would require considerable 
attenuation. Therefore the power it lieited to 5 Vatta for both the punp and 
Stokes latere. A good choice for thit later tay be a high repetition frequency 
NdtYAG later that has cw flash lamps and a Q- switched cavity. A typical later of 
thit type would have a pulse repetition frequency (prf) of 10 kHx and a pulse 
duration of 233 nanoseconds. In 10 Billiseconds 100 pulses would be obtained 
giving a dye later energy of 0.1 tl. The estimated CABS signal count would be 2 
x 10 3 photons for a 100 cn" 1 wide dye laser. 

The estinated signal magnitude depends critically on the laser power for 
nonlinear processes such as CARS. An optical fiber is expected to be capable of 
handling a higher intensity of pulsed radiation if the pulses are short enough in 
duration. The pulse duration must be shorter than the time required to traverse 
the optical fiber. For a 40 meter long fiber this tine would be about 200 
nanoseconds. Therefore, it is expected that slightly higher peak powers could be 
used in the case just described, but this gain would not be sufficient to offset 
the low signal. Nd:YAG lasers are commercially available which operate with muen 
shorter pulses (1C nanoseconds), using pulsed flashlamps; however the prf is 
restricted to 10-20 Hz. This low frequency is incompatible with the desired 
temporal resolution. 

Although the CARS approach would be more attractive with laaers intermediate 
between those two operating ranges, the CARS technique also has other disadvan- 
tages, aside from the requirement for through optical access. The nonlinear 
susceptibility depends on the Raman linewidth, which in the extreme conditions of 
the SS>£ pre-burner requires extrapolation of existing linewidth data. This 
would need to be verified in experiments. Another potential problem is frequency 
shifting, which for the cw approach may require retuning the Stokes laser to stay 
in resonance with the molecular transition. This can be accomplished through a 
variety of frequency stabilisation approaches, but makes the instrument more 
complicated. 
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Lilir Induced Fluorescence 


The laser fluorescence power ie « function of the laser intensity, the 
quenching rates of levels involved in the transition, the spontaneous emission 
rates and the frequency overlap with the absorbing molecular frequency (Refs. 3 
and 4). The fluorescence power varies linearly with laser power at low power, 
but reaches a constant value at high powers when the radiation is intense enough 
to saturate the molecular transitions. 

The expression for the fluorescence power, 8p, assume a sinple fora at 
saturation 


S F = hkp A 2| & V c N,j 


(3) 


where h is Planck's constant, Vp is the transition frequency, A 21 is the radia- 
tive transition rate, Q and V £ are the solid angle and volum respectively from 
which fluorescence is collected and Njj is the population of the lower 
(absorbing) state involved in the transition. For the pre-burner conditions the 
most populated lover level, J' ■ 4.5, in the notation adopted for spin-orbit 
coupling, has a population of 1 .0 x 10 7 ca~ 3 . The radiative lifetim is 0.803 
microsecond (Ref. 13). It is estimated that in 10 milliseconds, 1.1 x 10 4 
photons would be collected which compares closely with the signal for the 
equivalent H 2 spontaneous Raman system. 

It is not likely that the transition involved in the transition can be 
saturated because of practical limitations on the incident laser power. Lucht, 
Sweeney and Laurendeau (Ref. 14) using a dye laser pumped by s Nd:YAG laser 
observed fluorescence from the directly excited J' ■ 4.5, N' • 4 upper rotational 
level using a spectrometer tuned to the Pj(5) transition. At atmospheric 
pressure the power needed to scturate the transition was 150 kW (at - 3100A). 

This value can be scaled linearly to the pressure conditions of the SSME pre- 
burner. Therefore at the nominal conditions obtained in the pre-burner the 
saturation intensity would be about 55 GW corresponding to intensities in excess 
of 7 IW/cm 2 in optical fibers of 1 mm diameter or less. This intensity is far in 
excess of the levels considered for this application. It means that the photon 
count estimated from Eq. (3) is too high, the actual count attainable would be 
lower. 
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Technique Selection 

Within th« constraints imposed on optical accass to tht SS1C prc-burncr the 
selection of the optimum technique is clear. CARS and other coherent techniques 
afford only temperature measurement along a line in the pre-burner defined by the 
chord joining the too existing transducer ports. CARS and like techniques can 
yield large signal with unlimited laser power, but when the power is limited to 
levels amenable to transmission through an optical fiber, the signal level is no 
higher than spontaneous Raman scattering. Another potential problem of a lesser 
nature is associated with the Raman linewidth at these gas conditions. Analysis 
of CARS spectra would require a knowledge of the widths of individual transi- 
tions. Existing data do not apply to the pre-burner conditions; therefore addi- 
tional measurements of these would be required. 

Laser induced fluorescence signal levels are predicted to be far below the 
other techniques for equivalent diagnostic conditions, due to the low concentra- 
tions of fluorescing species, namely OH. The probing of molecular levels in OH 
would require a source emitting in the ultraviolet and more importantly, the 
transmission of uv radiation through an optical fiber. The absorption of fibers 
generally rises in the short wavelength region of the optical spectrum due to 
Rayleigh scattering, making the transmission of uv relatively poor. Since 
saturation of the fluorescent transition is not practical, a knowledge of 
quenching rates is also required. The use of a seed material could increase 
fluorescence levels, but may be difficult to implement. This is an unnecessary 
complication since the Raman technique appears capable of achieving the 
diagnostic design goals without seeding. 

The spontaneous Raman backseat taring approach alleviates all of the afore- 
mentioned objections. It appears practical to make temperature measurements 
using this technique. This has been indicated somewhat by the calculations shown 
thus far, but calculations shown in the next Section of this report carry the 
predictions further. One drawback, to be considered with Raman scattering is the 
relatively low collection efficiency imposed by the limited window size in the 
SS1C pre-burner. The very high density offsets this limitation however. 
Frequently in practical combustion measurements the background radiation due to 
flame emission limits the Raman signal-to-noise ratio. However the hydrogen 
flame exhibits extremely low visible emission in the laboratory. This is also 
expected to be the case in the pre-burner due to the fuel-rich condition. 

Detailed knowledge of the linewidth and line shift of transitions is not needed 
for Raman temperature measurements. One only needs to be satisfied that the 
transitions do not become so broad that there is appreciable overlap. This is 
not expected, as will be shown in the next section. The theory of Raman 
scattering in hydrogen is among the most well known, so that no theoretical 
hurdles are expected for the analysis of data. Finally and very significantly 
the Raman backseat ter ing approach appears to be amenable to implementation with 
optical fibers. Furthermore only one laser source is required, ameliorating the 
expense and practical problems associated with operating two lasers. 
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This section describes spontaneous Kmn scattering diagnostics, selected as 
the best technique for SSME pre-burner optical diagnostics. This section begins 
by describing the principles of Raman temperature measurements as applied to the 
fuel pre-burner. The theory of hydrogen Raman line intensities, and lineshapes 
is presented. The implications for diagnostics of the two possible nuclear spin 
orientations in hydrogen are identified and the results of calculations for the 
nominal 'power condition are presented. Following this the design of the optical 
diagnostic system is considered. The use of optical fibers is examined and the 
desired characteristics of other system elements such as the spectrometer, 
detector and laser are enumerated. Finally the perfotmance of the Raman 
diagnostic system is estimated, and evaluated in terms of program design goals. 


Principles of Raman Diagnostics 
Theory of Hydrogen Rastan Intensities 

The calculation of theoretical Raman spectra is straightforward. Raman 
spectroscopy is well presented in recent reviews (Refs. 5-7). In the usual 
prescription, the Raman scattered power Pg is related to the incident laser 
power by 


w Hi) 


(A) 


where do/dft is the Raman cross section, i^j is the population density, t is 
scattering length, and ft is the solid angle over which Raman scattered radiation 
is collected. This simple expression is inadequste for the expression of Raman 
scattering powers in hydrogen, because of the effects of the rotational-vibra- 
tional interaction. This factor which is most important for anharmonic 
oscillators, results from the perturbation of the vibrational energy levels by 
centrifugal distortion of the molecule in higher rotational levels. James and 
Klemperer (Ref. 15) first studied this factor, which is analogous to the Herman- 
Wallis correction to the transition dipole moments in infrared spectroscopy. The 
correction factors for H 2 and other atmospheric gases are described by 
Asawaroengchai and Rosenblatt (Ref. 16). These factors have been used in the 
calculations presented here. Recently Cheung et al. (Ref. 17) calculated more 
accurate centrifugal distortion factors. These result in a difference of less 
than 2Z for the tenth rotational level, therefore analytic expressions of James 
and Klemperer have been retained. 
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Raman scattering arises from the induced polarizability of the molecule, 
that it tha ralativa ratponta of tha alactron cloud and tha nuclai to tha 
electromagnetic fiald of tha incidant radiation. Tha polarisability ia a tanaor 
which can ba decomposed into an iaotropic part and an anisotropic part. For 
diatomic molecules tha anisotropic part is simply 0 • a ( - o^, whsrs is tha 

polarisability element parpandicular to tha auclaar axis and a ( is tha 
component parallal to tha axis. Tha isotropic part is «■ (a ( ♦ 2«.)/3. 

Tha intansity of a Raman molecular transition in terms of anargy flux 
seattarad froa unit voluaa par solid angla is given by: 


where I Q is the energy flux of the incident radiation, v is the frequency of 
tha scattered radiation, N y j is tha density of scattering aolaculas, v and J 
are the vibrational and rotational quantua numbers respectively and i t the 

square of the polarisability matrix element. The matrix element squared can be 
written in the form (Ref. 15) 


<a > * C b S (J) f (J) 


( 6 ) 


where C b is an intensity factor which depends on the rigid rotator polaris- 
ability, S(J) is the rotational line strength and f(J) is a correction for the 
centrifugal distortion caused by vibrational-rotational interaction, derived by 
James and Klemperer (Ref. 15). 

The line strength for J -> J ♦ 2 pure rotational Stokes (0-branch) 
scattering and Stokes vibrational-rotational S-branches is 


. ( J+ l) (J + 2) 


S(J, J+ 2) 


2 J + 3 


(7) 
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snd the line strength for J -> J - 2 transitions is (S-branch) 


S(J,J-2)* 


(J-l) J 
2J-I 


( 8 ) 


Joes and Klemperer showed that for pure rotational Stokes transitions 


f(J) 


$ I 
00 



(J*+ 3J + 3 )]* 


(9) 


where x* 8 e /r e 6 e *, ®e* ** the derivative of the polarisability anisotropy, 

r e is the equilibrium nuclear separation and B e and o> e are the rotational and 
vibrational constants respectively. For anti-Stokes transitions 


f (J) 


o 

00 




( 10 ) 


where here the notation has been changed somewhat in that J denotes the 
rotational quantum number of the initial state, so that the transition considered 
is J -> J - 2. 

For rotational-vibrational transitions the f correction factors are given 

by: 


S-branch 


[-ft)' 

4x( 

’it) 

i w e / 

(2J+3)] 2 

(11) 

0-branch 


[-ftr 

J{J + I) + Ax 1 

(B, 

) (2 J “ * >] * 

(12) 

Q-branch 


f(J>or * [ 1 * 

j+ 


(13) 
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The rigid rotor intensity factor it given in titaa of tha molecular 
poiaritability component a which art auanariaad in Tabla 111. 0 q ia tha 
polariaability aniaotropy for v • 0, r c ia tha equilibrium intarnuclaar radiua 
and 6 C * ia the firat derivative of the polariaability aniaotropy (d6/dr)r c . 


TABLE 111 

MOLECULAR POLARIZABILITY FACTOR, C fe (v - 0) 


Pure rotation: C^ ■ (7/30) 

Vibration-rotation C^ ■ (7/30) (r e 8 e ')^ (B e /o> e ) 


The correaponding relationahip for the Q-branch intenaity ia not expreaaible as a 
simple factor. James and Klemperer explicitly derived the expression for the 
case of Raman scattering being observed at 90 degrees to the direction of 
propagation of the exciting radiation, but the expression for backscatter can be 
written down from their equations. The corresponding expression for the Q-branch 
intensity in the backscatter direction is 


I®, =(!»k,)* (®l) f®, (j) n oj [(r, a',) 1 * rUHr e r') 8 /<2j + i)] (14) 


where 


/ tJ) - I J(J + I)(2J+I) _ 2 J+l + 4 J(J-H) JtU)+ JK* I) 

15 (2 J+3)(2 J * I) 9 15 2 J+l 15 


(15) 


and 




J(4J 2 + I) 


(2J + I)(2J-I) 


( 16 ) 
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The parameters used to evaluate those sxprsssions art summarised in 
Table IV. The centrifugal distortion factors have been evaluated for the 
rotational energy' levels populated in the SSME pre-burner flow conditions. 
Rotational levels up to J ■ 5 or 6 are significantly populated. The centrifugal 
distortion factor is 1.243 for an 8-branch transition from this level and 1.129 
for an Q-branch transition. These contributions will be important for accurate 
diagnostic measurements. Rosenblatt and co-workers (Ref. 17) made more accurate 
ab initio calculations of these factors and found differences with Eqs. (9-13) 
that are important for rotational levels J > 8, and are therefore significant 
only in hotter flames T > 1800. 


TABLE IV 

POLARIZABILITY MATRIX ELEMENT DATA 


Quantity 


Value 


Reference 


*0 

0.314 x 10" 2 ** cm 3 

16 

X 

0.38 

16 

a' 

1.230 x 10" 16 cm 2 

18 

Y* 

1.02 x 10" 16 cm 2 

16 


0.741 x 10" 8 cm 2 

16 


Hydrogen Lineshapes 

The linewidths and frequencies of Raman lines are affected by atomic and 
molecular collisions. An increase in pressure can cause broadening and shifting 
of the Raman lines of hydrogen. Considering the extreme pressures existing in 
the SSME pre-burner it is important to determine the magnitude and impact of 
these effects on optical diagnostics. 
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Pressure broadening and ahifting has bean tha subject of numerous 
investigations because intermolecular forces can be atudied by tonsuring the 
effects of pressure and temperature on spectral linewidths (Ref. 19). At low 
pressures (less than 0.1 atmospheres) and moderate temperatures, the Reman 
linewidth is Doppler broadened and its shape is Gaussian. At pressures greater 
than 10 atmospheres the Raman linewidth is pressure broadened, the linewidth 
increases linearly with gas density, and the lineshape is Lorentsian. The 
linewidth is proportional to Bp, where B is the optical (Raman) broadening 
coefficient and p is the density, normally in amagats. In the range of densities 
between the Doppler limit and the pressure broadening regime, line narrowing, 
named for Dicke who first described it (Ref. 20) occurs in hydrogen. Doppler 
broadening for a scattering event arises from the recoil momentum imparted to the 
molecule by the scattering photon. If an elastic collision occurs simultaneously 
with the scattering event, some of the momentum may be taken up by the colliding 
molecule; and the scattered photon will have a smaller Doppler shift than if no 
collision had occurred. Pressure broadening becomes dominant when the collision 
rate is so large that a substantial number of inelastic collisions occur, which 
change the internal state of the molecule. Dicke narrowing takes place when the 
mean free path for elastic collisions is comparable to the wavelength of the 
scattered photon and is less than the mean free path for line broadening 
(Ref. 21). 

As the density is increased another narrowing mechanism is manifested in the 
Raman lines of most molecules. Collision or motional narrowing occurs when 
closely spaced Raman lines broaden into one another, and the transitions can no 
longer be considered to be individual transitions. This occurs at quite low 
pressures in common atmospheric gases. Due to the very large frequency spacing 
of lines in hydrogen it does not occur until densities corresponding to several 
thousands of atmospheres have been reached (Ref. 22). 

The Raman linewidths corresponding to rotational and vibrational transitions 
have somewhat different broadening coefficients dee to differences in the effec- 
tive collision mechanisms which affect the two types of transitions. This has 
been discussed by Srivistava and Zaidi (Ref. 19). The reason for the difference 
is that rotational linewidths are the result of elastic and inelastic collisions. 
Therefore the rotational linewidth is broader than the vibrational linewidth. 
Table V shows the broadening coefficients at 300 K. There is very little data on 
the effect of temperature on the Raman broadening coefficients. Lallemand and 
Simova (Ref. 23) measured the broadening coefficient from 300 to 425 K. Their 
data indicate a T 1 * 25 temperature dependence for the Q-branch transition Q(l). 
Assuming the broadening coefficient varies with this temperature dependence and 
taking an average broadening coefficient at 300 K of .0025 cm~ 1 /amagat the 
estimated linewidth at SSME pre-burner conditions (5500 psia, 1000 K) is approxi- 
mately 0.5 cm -1 . This is less than the expected spectral resolution of the 
spectrograph to be used and therefore the actual linewidth is not critical. 
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RAMAN BROADENING COEFFICIENTS 


Rotational 


Lincwiith x 10" 

3 

Transition 

(Ref. 24) 

(Ref. 25) 

(Ref. 26) 

8 0 (0) 

2.80 

3.1 

3.0 

8 0 (D 

3.47 

3.8 

4.6 

S 0 <2> 

? .47 



8 0 (3) 

2.53 



S 0 (*) 

1.53 



Vibrational 


Linewidth x 10" 

3 

Transition 

(Ref. 27) 

(Ref. 28) 

(Ref. 26) 

Qoi«» 

2.32 




1.40 

1.68 

1.90 

Qoi<2> 

2.53 



9oi<3> 

3.66 




Similarly the effect of preeture on the frequency of Reman transitions has 
been measured in hydrogen (Ref. 19). The shift is approximately -0.2 cm -1 /100 
amagats at 300 K (Ref. 23) for the Q(l) transition. Lallemand and Simova 
(Ref. 23) found the linear pressure shift coefficient to vary linearly with temp- 
erature from 80 to 500 K, becoming positive for temperatures above 420 K. 
Extrapolating these data to the operating conditions of the SSME pre-burner, the 
shift is very roughly estimated to be on the order of +1 cm" 1 . Rotational shifts 
have been found to be smaller (Ref. 26). There is little impact on optical 
diagnostics expected due to these shifts. 

The Effects of Ortho and Par a Hydrogen Modifications 

Hydrogen has a nuclear spin of 1/2. Different orientations of the two 
nuclear spins with respect to one mother give rise to two molecular modifica- 
tions designated as ortho and para. In the ortho form, the nuclear spina are 
aligned parallel, and in the para form the spins are anti-parallel. The states 
with even rotational quantum number J belong to the para-modification, while the 
odd-J states belong to the ortho-modification. The relative equilibrium composi- 
tion depends on the temperature. At room temperature hydrogen is 75 percent 
ortho-hydrogen, but at the normal boiling point, 20.4 K, hydrogen is 99.79 
percent para-hydrogen. 
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The conversion of ortho to pars is slow in the sbssnet of catalysts at lov 
to moderate temperatures so that it is quits possible to have a composition that 
differs greatly from equilibrium. The conversion of hydrogen from the orth-'-form 
to the para-form is accompanied by the release of 1.06 kJ/mole of heat at the 
normal boiling point (Ref. 29). The heet gs terated in this reaction is compar- 
able to the heat of vaporisation, so that this reaction can account for boil-off 
losses in liquid normal-hydrogen (75 percent ortho-hydrogen) of about IX per hour 
(tef. 29), which ere unaccepteble in most applications. The boil-off losses in 
the pure para-modification, however due to the absence of conversion, can be 
minimised by good tank design practices to levels below 1 percent per dsy. 
Therefore hydrogen is almost always converted catalyticelly to the para-modifica- 
tion during liquefaction. 

The direct spontaneous transition from ortho to pars through spin-inter- 
change is strictly prohibited (Ref. 30). 8elf-conversion takes place through 
atom-exchange: 


H 4 p-H 2 — o-H 2 4 H 


(17) 


This reaction does not violate the prohibition of spin interchange because a new 
molecule is formed in the reaction. The atoms are supplied in this reaction by 
dissociation of molecular hydrogen. Spontaneous conversion was studied by Farkas 
and Farkas (Ref. 31). Subsequent measurements suggested that the earlier results 
were affected by small quantities of oxygen in the hydrogen stream (Ref. 32). 

The reaction time t, is related to the reaction rate constant k (Refs. 33, 34), 
and the equilibrium constant K, by the relation 


t * 


KK [hj] 


(18) 


The time calculated from this relation at the conditions existing in the SSME 
pre-burner is approximately 2 sec. This assumes that the mechanism which has 
been confirmed at low pressures does not change over several orders of magnitude 
in the pressure. Oxygen present in the SSME pre-burner will affect the H-etom 
concentration and also catalyse the converion of para-hydrogen into ortho- 
hydrogen (Ref. 32), so that actual conversion rates will be fester than the self 
conversion rate calculated here. 
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The characteristic reaction time for self conversion is longer then the 
residence tins in the pre-burner. Therefore it is possible that the hydrogen 
ortho-para composition is not in equilibrium. This mill effect the predictions 
•of the rotational state populations in the measurement volume. 8ince the pare 
modification dominates at low temperature it means that the even rotational 
levels will be more highly populated than predicted from equilibrium, tending to 
minimise the differences in the two manifolds. The impact of these considera- 
tions is on the analysis of the date. When the relative populations are 
analysed, the even and odd rotational levels must be considered separately. Both 
rotational manifolds should yield a valid temperature, and this ran serve as a 
check on the consistency. The total populations will not be greatly affected so 
that the predictions of Raman scattering signals still remain a valid estimate of 
the system performance. 

Results of Raman Calculations 


Calculations have been performed based on the theoretical development des- 
cribed for the effects of anharmonicity on the hydrogen Raman intensities. The 
calculations are illustrated in Figs. 2 and 3. The first of the figures depicts 
the calculated Raman spectrum of hydrogen at 38 )fa (5500 psi) and 1000 K corres- 
ponding to the nominal conditions in the fuel pre-burner. The Raman intensities 
are shown as a function of Raman shift from the exciting wavelength, 4880 A for 
the calculation shown. The pure rotational Raman transitions corresponding to 
the 0 (AJ ■ -2) and S (AJ - +2) branches are shown with negative and positive 
frequency shifts respectively. The vibrational Q-branch (Av ■ 1, AJ ■ 0) is 
shown near 4160 cm' 1 shift, with associated vibrational 0 and S branches near it. 
The various transitions arise from individual rotation transitions. The rota- 
tional S branch is nearly 2 times stronger than the vibrational Q branch at these 
conditions, indicating that rotational Raman would produce better system perfor- 
mance. 

The sensitivity of the distribution of these rotational lines to gas 
temperature is shown in Fig. 3. The rotational quantum number of the initial 
level of the transition is shown in the bottom spectrum in the figure. The 
sensitivity of the distribution of the lines is shown. Bow temperature is 
deduced from Reman intensities is described in the next section. 

Appendix A contains the output from a computer program based on the formula- 
tion described earlier. This tabulation gives for the verious transitions the 
Raman shift, the wavelength of the shifted line, the molecular populations and 
the Rasun intensities. 
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Um p Temperature Measurements 

The sensitivity of the hydrogen Kasrn intensities to temperature is shown in 
Fig. 4. In practice, the gas temperature is determined from measured Raman 
intensities by appropriately normalising and plotting them on a semi-logarithmic 
plot as a function of the rotational energy of the lower level involved in the 
transition. The principle of this procedure is illustrated in Fig. 4, which 
shows the normalised Raman intensity distribution for four temperatures from 600 
to 1200 K. normally the Raman aignal would be in units of photon or photoeiee- 
tron counts from a detector. These signal intensities are normalised by dividing 
by v 3 g(J)S(J)f(J)/R(v) where R(v) accounts for the spectral response of the 
optical system, a factor which must be determined empirically with a standard 
spectral lamp. The temperature is determined from the slope of these lines, m, 
since m - hc/kT. As a practical matter the line intensities may need to be 
determined by integrating across the spectral width of the lines on the multi- 
channel optical detector. 


System Design Considerations 

The system for the SSME pre-burner optical diagnostics is conceived as 
composed of a high power laser, an optical head mounted directly on the pre- 
burner, and a spectrograph/optical multichannel detector (OMD) combination. The 
laser and the spectrograph/OMD would be located remotely, and would be joined to 
the optical head by optical fiber links. The use of optical fibers is crucial to 
this application because of the severe environment surrounding the operating 
engine. The theoretical limits to propagation of high power laser radiation are 
discussed in the next section. Following that is a discussion of the selection 
of the laser source. Lastly the requirements for the spectrograph/detector 
combination are evaluated. The elements of the diagnostic system are shown in 
Fig. 5. 

Power Handling Capability of Optical Fibers 

The power handling capability of low-loss optical fibers is generally 
limited by nonlinear processes occuring in the fiber for visible and near 
infrared radiation. The surface quality and finish of the fiber end also can 
affect the critical power for fibers carrying ultraviolet radiation. 

The theoretical power handling capacity of small diameter (single mode) low- 
loss optical fibers was investigated by Smith (Ref. 35). He found that stimu- 
lated Raman and stimulated Brillouin inelastic scattering determine the practical 
power capacity of fibers at high power densities. In Raman scattering, the pump 
laser radiation field interacts inelastically with the nuclear vibrational mode 
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to product t frtqutncy which it down-shifted by tht nuclttr vibrational 
frequency. In stimulated Raman scattering, tht optical wave itatlf producta a 
cohtrtnt wavt in tht nuclttr vibrational coord inatt which intaracta with tht 
tranaaitttd optical boat. Stiwulattd Brillouin ariata from tht cohtrtnt teouatic 
wavt drivtn by tht optical fit Id. Tht difftrtnct batwaan Raman scattering and 
Brillouin acattaring it that tht former excites an optical phonon while tht 
latter axcitaa an acouatic phonon. For forward acattered Raman, the principal 
affect ia to cauat the pump frtqutncy to be thiftad and broadtnad* Backward 
processes, either Raman or Brillouin) can raault in atvart attenuation of the 
forward traveling optical wavt, due to tht stimulated backward wavt, and there** 
fort poaa a wort atvart linitation to fiber power tranawiaaion. 

Tht criterion that Smith chose to determine the critical power for nonlinear 
proctasts waa taken arbitrarily to be that power for which the optical field 
generated by the nonlinear interaction ia comparable to the incident field. 

These optical processes depend exponentially on the pump power level, ao that a 
reduction of the punp power level of 1 dB, reduces the etiaulated scattering 
level by 20 dB (Ref. 33). Therefore the calculated critical power level serves 
aa a useful benchmark. 

Saith aaauaed that the prop laser radiation propagates in the fiber aa 
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p p (2) * Pp (0) exp ( - a p z) 


(19) 


where Pp(t) ia the punp power at position a along the fiber, Pp(0) is the 
incident power and Op is the attenuation coefficient at the pixap wavelength. 

This means in effect that the puap field is not affected by the nonlinear inter- 
action. Based on the above definition of critical power, Saith derived the 
following expression assuming the fiber length ia much greater than the absorp- 
tion length 



2 




( 20 ) 


where h is Planck's constant, v is the Stokes shifted frequency (Raman or 
Brillouin), y is the nonlinear gain coefficient, A ia the fiber cross-sectional 
area, and Ak is the Raman or Brillouin linevidth (FWHM). The factor f and the 
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exponent n depend on the etimileted process. For forward Raman scattering 
n ■ 3/2 end for backward nonlinear processes n ■ 5/2. The factor f is 1 for 
Raman scattering (forward or backward) and for Brillouin scattering f * kT/hv a 
where v a is the frequency of the acoustic phonon. 

In order to evaluate the critical power it is necessary to know the gain and 
attenuation coefficients for the optical fiber at the wavelength of interest. 

For wavelengths in the visible away froa the hydroxyl absorption bands near 0.7- 
1.4 micron, the attenuation coefficient can be estimated from Rayleigh scattering 
(Ref. 36). The estimated abeorption coefficient is 

a p = 2.6 x 10 '• \p' 4 (21) 


where the wavelength X p is in microns and the absorption coefficient is in 
units of cm -1 . The absorption coefficient in this formula corresponds to an 
attenuation of 1.7 dB/km at 0.9 microns. Stolen (Ref. 37) has found that the 
Raman gain coefficient scales linearly with pump frequency while the Brillioun 
gain coefficient is almost independent of pump frequency. Taking the gain 
coefficients from Smith (Ref. 35) then the appropriate gain coefficients are 


G„ = 5X 10 ' ,0 X P * 4 cm/w (22) 


G B = 3 X 10** Cm/w (23) 

where X p is the pump wavelength is microns (Ref. 37). Smith found that the 
solutions of Eqs. (20) are not sensitive to the Stokes frequency, the linevidth 
or the factor f. The solutions yield (Ref. 35) 


Pc s 


16 AQp 
g r 


watts 


(24) 
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P c « Afl f Watt* (25) 

6 i 

The critical intensity then it given by 


Xj 8 8.32 x to 4 X P ‘ 8 W/cm* 

I® * 1.82 x I0 4 X P ' 4 W/cm 2 


(26) 

(27) 


where again the wavelength is in microns. These foraulae are derived assuming 
that the fiber is longer than the absorption length, i.e. 1/oip. For radiation 
at 0.4 microns wavelength, this length must be in excess of 100 meters. This 
length is much longer than the length contemplated for SSME prr • „ jrner optical 
diagnostics. In this case the formulae for the critical powers are modified by 
replacing 1/% by the effective length, *eff : 


eff 


i- exp (-a p I ) 


(28) 


and 1 is the actual length of the fiber. The critical powers for stimulated 
Raman and Brillouin scattering are shown together in Fig. 6. A strong dependence 
on wavelength is predicted. This is due to the sharp increase in absorption for 
shorter wavelengths which attenuates the stimulated wave increasing the power 
that the fiber can handle before these processes become severe. For wavelengths 
as short as 0.4 microns the transmission for the 40 m length considered is 
greater than 80 percent. 

The analysis above is relevant when there is no reflection from one of the 
ends of the fiber. Labudde, Anliker and Weber (Ref. 38) pointed out that when 
both reflectivities are finite an optical resonator is formed with a threshold 
pump power P ( given by 


/9 P t L «ff \ , . . 

r,r 2 exp \ 1 exp(-aL) 8 1 


( 29 ) 
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ESTIMATED CRITICAL POWER IN FU8ED SILICA FOR STIMULATED PROCESSES 



WAVELENGTH, nm 


• 3 — 6 — 129—3 
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Labudde at al. alto pointad out that whan a laaar ia coup l ad directly into tha 
and of tha fibar tha laaar output mirror contributes to tha affactiva reflectiv- 
ity of tha fibar. In axpariaanta thaaa authors usad a Clan-Thoapaon priaa and 
quartar wava plata combination to dacoupla tha laaar froa tha fiber, noting that 
this ia only partially affactiva bacauaa tha backscattarad light, which is 
unpolariaad dua to polarisation acraabling in tha fibar, is attanuatad by only 50 
pareant in the isolator. The critical threshold intensity for Raaan and 
Brillouin scattering calculated froa Cq. (29) is shown in Fig. 6 for two and 
reflectivity combinations, by curves labelled "Brillouin oscillation threshold". 
Tha lower curve corresponds to tha case whan there is no coupling isolation and 
shows its importance. The decrease in critical intensity with wavelength is not 
as severe as before because of the dependence on the affactiva length. Tha 
critical Raaan intensity increases at longer wavelength because of tha lower gain 
coefficient. Since the Brillouin gain is independent of wavelength, its critical 
intensity just decreases with wavelength. Comparing the curves shown in Fig. 6, 
it is seen that the critical intensity for wavelengths around 0.5 aicrons is 
limited by the feedback from the ends of the fiber by Brillouin scattering. When 
the effective fiber end reflectivity is larger than 0.30, the critical intensity 
is even lower. 

These considerations apply to the transmission of a single transverse mode 
through the fiber. The power liaitation for a single mode fiber is on the order 
of tens of milliwatts. This is too small for SS1C pre-burner optical diagnos- 
tics, therefore the use of aultimode fibers is envisioned. The critical power is 
correspondingly higher for multimode fibers (Ref. 35). Table VI shows the 
critical power per transverse mode calculated for an intensity limit of 40 
kW/cm 2 , which Fig. 6 indicates is reasonable for Brillouin scattering. The 
values shown are per transvese mode. Typically an argon ion laser such as 
envisioned for these measurements would excite several transverse Bodes, so that 
the use of total power 5-10 times these figures seems reasonable. 


TABLE VI 

CALCULATED CRITICAL POWER PER TRANSVERSE MODE 


Diameter 

P 

(micron) 

JSi 

50 

0.784 

100 

3.14 

150 

7.08 

200 

12.6 
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Smith point* out that th* critical power for fibara can ba highar for pulsed 
laaar aourcaa. The analyai* abova i* applicable only to ew laaer aourcaa. When 
the pulae duration ia leaa than tha tiaa required to travaraa the fiber than th* 
laaar pulaa doe* not interact with tha generated stimulated radiation along the 
full length of tha fiber. Tor tha length of fiber anviaagad for thi* application 
th* pula* duration would need to be laaa than about 100 nanoaaconda for thia 
affect to be important. Thi* doe* not lead to large incraaaaa in tha critical 
power for pulaad laaar aourcaa though. It ia not lihaly to wake pulaed laaar 
aourcaa wore attractive becauae th* highar intanaitiaa they afford are atill not 
practical in fibara. 

Recently Brunei at al. (Ref. 39) demon* t rated the transmission of argon ion 
laaar radiation, like that being conaidered for thia application, up to 18 Watt* 
through a fiber a* email aa 150 micron* with meaaured power tranamiaaion of 
approximately 60 percent. These level* ahould be adequate for SSME Raman back- 
acatter diagnostics. Teat* at UTRC confirmed the measurement* for a low power, 
0.7 W, argon laser operating on the 4880 A line. These measurements support the 
theoretical astimatea described above. Further tests would be required to 
confirm the limitations of particular fibers. 

Spectrograph/detector 

The use of the optical fiber in the transportation of the collected Raman 
signal requires some consideration in the design of the spectrograph. The 
spontaneous Raman radiation collected by the optical head, and focused into the 
optical fiber, being incoherent, would emanate from the fiber at the characteris- 
tic numerical aperture of the fiber. Typical fibers have a numerical aperture of 
0.2 to 0.4 corresponding to f-numbers from 2.5 to 1.25. This is faster than most 
fast general purpose spectrographs which are typically f/8. Recently however 
some spectrographs have been designed employing concave gratings which are f/4. 
The fiber being envisioned for this application is a 50 micron core graded-index 
fiber with a numerical aperture of 0.14. This nominally requires a f/3.5 spec- 
trometer. A spectrograph has been designed for use in CARS tests at UTRC which 
is f/3.5, employing a 20 cm diameter concave holographic grating with a groove 
spacing corresponding to 2000 grooves per mm. Thia particular instrument with 
this grating has a dispersion of 5 A/ am. 

The required spectral dispersion depends on the Raman transitions being 
observed. A larger spectral working range is needed for rotational Raman 
scattering than for vibrational Q-branch diagnostica. The scattered wavelengths 
corresponding to rotational transitions span 285 Angstroms, from 4965 to 5250 A 
for 4880 A laser excitation for the SSME gas conditions. In order to disperse 
this spectrum across an optical multichannel detector face of 25 mm requires a 
dispersion of 11 A/mm. However one may choose to disperse over half the active 
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tru to ovoid edge effects so that a diaparaion of 25 A/mm way b a praftrrad. 

Thia diaparaion can be achieved by proper aelection of the grating groove 
apacing. The vibrational hydrogen Q-branch spana a wavelength range from 6123 A 
to 6063 A for 4880 A excitation. Thia would require a diaparaion of 5 A/m to 
diaperae the apectrum acroaa half the detector face. 

For unit magnification in the spectrometer, a 50 micron fiber coupled 
directly into the entrance would limit the reaolution to two 25 micron wide 
pixela on the detector. Since the dieperaion optimised for rotational Raman 
correaponda to 0.1 cm"* per micron , the inatrumental apectral reaolution would be 
5 cm"*. With aome defocuaing, which ia characteriatic of diode arrays, this 
reaolution may be degraded' to 8 cm" 1 . The apectral reaolution for vibrational 
Raman diagnoatica would be higher. 

The detector enviaaged for thia application ia a aelf scanning diode array. 
Thia detector haa good optical response at the wavelength of use, corresponding 
to 0.2 counts per photon. Each detector element or pixel can be individually 
accessed, processed and converted in 16 microseconds. Unwanted pixels can be 
bypassed at 0.5 microseconds/pixel. Therefore a spectrum consisting of 6 or 7 
lines can be scanned in 2 milliseconds when allowance ia made for reading 10 
pixels around the Raman lines. The diode array is much leas sensitive to vibra- 
tion than a vidicon detector, and does not suffer from blooming. The detector 
does have a limited dynamic range requiring care not to saturate the detector. 
This limits the dynamic range for one scan to 2000 counts. 

Figure 7 depicts the design of the spectrograph. Thia is similar to the 
spectrograph used for CARS measurements. A fiber transmits the collected radia- 
tion to the entrance of the spectrometer where a prism may be used to afford 
access to the grating entrance pupil. The light from the turning prism diverges 
to fill the concave grating which refocuses and diperses it on the face of the 
diode array. The grating is designed to work slightly off axis as shown, and 
produce a flat field. Minimum grating tuning is generally required, dispensing 
with the need for a sine bar and precision lead screw. The Raman spectrum 
dispersed on the diode is illustrated in the top of the figure. 

Laser Selection 


The optimum laser for Raman diagnostic is determined by the optical proper- 
ties of the fiber, the response of the detector at the Raman shifted wavelength 
and the characteristics of Raman scattering. Raman scattering, as can be seen 
from Eq. (5), depends on the Stokes frequency to the fourth power, favoring the 
use of short wavelength lasers. The response of the detector is relatively 
insensitive in the range of 5000 to 8000 Angstroms falling substantially for 
wavelengths outside this range. The fiber critical power for the feedback 


FAST RAMAN SPECTROGRAPH 




MS-956181F 


Halted case increases for wave lengths leea than about 4500 A, however the 
fraction of the radiation transmitted decreases due to the fiber attenuation. 

The argon ion laser seeas to be a good choice, with the 4880 line being favored 
over the 5145 A line. 

Optical Head Design 

The optical head is conceived to be bolted to the 8 SIC pre-burner. It would 
serve to focus light froa the laser into the coabustor and to collect the Raman 
backscattered radiation and send it to the spectroaeter/detector. The optical 
head design concept is shown in Fig. 8. Beyond the window in the coabustor are 
the actual elements of the optical head represented scheaatically. First is a 
device to slew the line of aeasureaent across a segaent in the chaaber. This 
device is shown represented by a series of optical wedges aounted on a wheel 
which successively indexes wedges of different angular deviation into the optical 
path. The wedge wheel is rotated by a stepper motor. The next optical eleaent 
is the combined focusing and collecting lens. This is represented as a single 
eleaent but nay contain a set of optical eleaents. Following this lens is a 
thick glass filter which is shown also holding a prisa used to direct laser 
radiation to the focusing lens. The glass filter is used to block unshifted 
laser radiation froa the collection systen. A central obscuration disk slso nay 
be incorporated into this filter. Individual lenses are shown collimating the 
incident laser radiation and focusing the Raaan radiation into the optical fiber. 
Finally shown are holders which provide adjustment of the optical fibers for 
alignment. 

The possible angles that the aeasureaent axis can be slewed through are 
United by optical access into the coabustor. A relatively thick window would be 
required to contain the pressure. This sets the distance that the slewing device 
must be back from the pre-burner. Off-axis scanning nay be restricted by this to 
approxiaately 5 degrees, before vignetting at the window severely decreases the 
Ramsn signal. 

The central obscuration disk is used to determine the spatial resolution of 
the coaxial system (Ref. 40). The actual resolution for a given obscuration 
diameter is best determined empirically. The effect of the blockage is not 
expected to be too severe because of the large area of the off-axis annulus. 

The laser beam is shown as a small diaaeter pencil of light inserted into 
the center of the optics by a small turning prism. The laser beam can be focused 
into the fiber by a relatively long lens, resulting in a small solid angle 
entering the fiber. This small solid angle is preserved if the fiber undergoes 
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very gradual bands in transit to tha optical head so tha radiation travels 
through the fiber in low order modes. If there many taicrobends the low order 
■odes couple to higher order mode which diffract wore on exiting the fiber. 


Raaan Backscattering Systea Estimated Performance 
Eagan Signal Levels 

The performance of Raman temperature diagnostic system can be measured in 
terms of the performance criteria or design objectives. These specify a spatial 
resolution of 0.5 to 1 cm, time response of 10 milliseconds, and a temperature 
measurement accuracy of 10 R, or 6 K. The performance of a Raman system can be 
estimated from the calculations of the normalised Raman intensities and the 
estimates of laser power that can be transmitted through the optical fiber. The 
estimates are also subject to the constraints imposed by the optical access in 
the SSME pre-burner. The transducer port is 1.1 cm in diameter and the combus- 
tion chamber is about 25 cm in diameter. Taking an optical aperture of 1 cm and 
a depth of 20 cm, the solid angle over which Raman backscatter can be collected 
is .002 steradians. The number of photoelectrons counted in the optical multi- 
channel detector is given by 


N 


e 

R = 


p i* * c c 

hck s 


(30) 


where Pg is the incident laser power, t is the spatial extent of collection, 
i.e. the resolution, 1^ is the normalized Raman intensity per cm per steradian, 
t is the temporal resolution, c c is the collection efficiency, rip is the 
quantum efficiency of the optical detector and k, the Stokes photon frequen- 
cy in wavenumbers. The factors in the denominator give the energy per Stokes 
photon. This equation is used to estimate the nuc' er of photoelectrons counted 
due to a single rotational Raman transition. The collection efficiency is the 
resultant of all Fresnel reflection losses and transmission losses in the Raman 
system. The transmission for the fiber would be about 85 percent for a fiber A0 
m long and having an attenuation of 20 dB/km at 5000 A. For temperatures in the 
range of 1000 to 1600 K the calculated normalised Raman intensity is 2 x 10~ 9 to 
2.3 x 10" 10 cnT l sR~ l for transitions involving the lowest six rotational levels 
which are moat populated. Taking the collection efficiency and the quantum 
efficiency each to be 0.2 then Eq. (30) can be expressed as 


Nr - 6x»0 4 



P l M 


(31) 
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With an incident laser power of 4 Watts, and a spatial resolution, 1 • 1 cm, then 
the number of OK) counts is 2.4 x 10* 5 t. 

The diode array was estiaated to be capable of being scanned in 4.5 as. In 
this period the signal would be about 1000 counts. This is within the linear 
range of the detector. The acquired counts froa previous diode scans can be 
added to effectively average over a longer period of tiae, decreasing the 
teaporal resolution. 

Teaperature Resolution 


An estiaate of the nuaber of signal counts required to achieve the teapera- 
ture resolution objective can be aade froa the expression for the rotational 
populations. It can be shown that the teaperature resolution is given by 



where it the population of the level corresponding to the quantum numbers 
v and J. The photon count is directly proportional to the populations lyj, 
therefore the populations in this expression can be taken to be the signal count 
as well. A 6*K teaperature resolution requires a signal-to-noise ratio of 120. 
For shot noise liaited statistics then the required nuaber of photon counts is 
15000, which according to the estiaate above would require 63 asec, or somewhat 
longer than the design 10 asec. This estimate is based on the least populated 
rotational level. The highest populated level can be nearly an order of magni- 
tude higher, for which the teaporal resolution criterion would be satisfied. In 
actual analysis, the temperature would be obtained from the distribution of Raman 
intensity, therefore the teaporal resolution would be expected to fall in between 
these two liaits. If tests indicate that more power can be propagated through 
the fiber then the teaperature resolution can be correspondingly higher. 

Potential Background Interferences 

These estimates give the expected signal level in absolute terms. It is 
also important to examine the Raman system performance with respect to background 
interferences. No optical measurements have been done in the SS)£ pre-burner, 
therefore the radiation background can only be estimated. 

Hydrogen and oxygen burn with a flame of very low visible luminosity. 
Therefore the most important background interference to Raman measurements is 
likely to be thermal radiation from either particles in the flow, such as from 
erosion, or from hot surfaces adjacent to the flow. In normal operation both of 
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these sources would be expected to be quite low. To determine the reletive 
magnitude of thermal background for conditions likely to obtain in the SSME pre- 
burner estimates of the possible perticle loadings have been made. The model is 
a aono-disperse distribution of particles of disaster d, with e nuaber density of 
n. 


Hottel and Sarofia (Ref. 41) derived an expression for the theraal radiation 
froa a collection of particles, by relating the eaissivity to the ebsorptivity. 
They showed that the eaissivity c, is given by 


€ * |- exp(-COi) 


(33) 


where c is the concentration of particles, a is the absorption area of the 
particles and 1 is the optical path length. For a nuaber density of 10 4 
particles per cm 3 , having a aean diameter of 10 aicrons and a path length of 25 
cm, the eaissivity is calculated to be 0.28. For particles at a temperature of 
1000 K the calculated grey body emission intensity is 5 x 10“ 11 watts/cm -1 cm -2 
sR" 1 . It should be noted that this particle loading corresponds to a high rate 
of erosion (70 gm/sec for a metal of density 4.5 gm/cm 3 such as titanium), and it 
is quite likely that this le\ .1 of erosion would be detected in aoae other manner 
and the test terminated. 

The thermal radiation intensity is a strong function of particle tempera- 
ture, whereas Raman scattering is much less sensitive to gas temperature. This 
is shown in Fig. 9 which shows the calculated spectral irradiance for Raman 
scattering with 10 4 Watt/cm 2 incident laser intensity and the calculated thermal 
emission for two emissivities as a function of temperature. The Raman irradiance 
is calculated to exceed the thermal radiation for most particle loading at 
temperatures below 1500*K, therefore the theraal radiation from particles is not 
expected under normal operation to be an impediment to Raman measurements. 

The pre-burner walls would have a larger eaissivity than a collection of 
particles, but the temperature is expected to be considerably lower. The strong 
temperature coefficient implies that the radiation from the wall would also be 
expected to be low compared to the Raman scattered radiation. 
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CONCLUSIONS 


The feasibility of rwot« optical diagnostics in tha apaca shuttle vain 
angina fual pr a- burner has baan invastigatad. Diagnostics baaad on Raman back- 
acattaring appaar to ba tha auparior technique for tha aadlua, offaring advan- 
tages over fluorascanca and coharant nonlinear techniques, such as CARS. No seed 
or additive apeciea and no Modifications to tha pra-bumer are needed. Raman 
diagnostics appaar to ba aaMnable to implementation with optical fibers, which is 
required by tha constraints of operating in tha environment of a high performance 
rocket angina. 

Tha use of Raman backscattering, based on theoretical predictions, appears 
to have a reasonable chance of meeting diagnostic system design goals. Again 
based on calculations, it appears feasible to transmit 5 to 10 Watts of argon 
laser radiation through a multimode optical fiber 50 to 100 microns in diameter 
to excite the Raman radiation. This level is predicted to be sufficient to 
achieve the program design goals of 6 K temperature resolution, 1 cm spatial 
resolution and 10 milliseconds temporal resolution. Some tradeoffs may be 
required when the effect of obscuration, needed to achieve the requisite spatial 
resolution in the coaxial geometry, is taken into account. This effect is 
difficult to quantify theoretically in simple terms. Obscuration also may 
sharply affect off-axis temperature measurements, requiring relaxation of the 
design goals. The difficulties normally associated with Raman scattering: weak 

signal strength and interference due to background radiation are not expected to 
be problematic due to the very high density in this application, and the low 
flame luminosity expected in the fuel-rich hydrogen- oxygen flame. There appear 
to be no theoretical hurdles or difficulties to analysing Raman spectra. 

Potential problems which are difficult to quantify in a paper study pertain 
to the optical quality of the high pressure combustion medium. There are no 
measurements of the optical quality of the flow. Strong density gradients due to 
mixing, combustion pressure waves, and two phase flow in the propellant spray 
could pose severe problems. In the worst case, these effects could make it 
impossible to determine the precise point of measurement, or if severe enough, 
prevent collection of the backscattered radiation. On the other hand, the 
dominant species in the flow, hydrogen, has a very low Gladstone-Dale constant so 
that at this pressure this medium is more suitable than others. In addition 
Raman scattering, being an incoherent technique, should suffer least from medium 
optical problems. Hot gas tests are need to assess the magnitude of these 
difficulties. The design of the window is expected to present a challenge. The 
window must provide a safe, positive seal of the combustion chamber, at the same 
time withstanding thermal shock, and the radiation heat load from the combustor. 
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APPENDIX A 


Th« calculated Raman fraquanciaa and intansit&aa art tabulated in thia 
appendix. The aolecular parameters give the energy level constants, apin statis- 
tical veightt t aolecular polaritability and incident laaar frequency. Theae 
constants can be changed to reflect a new set of energy constants! or updated 
polarisabilities. The term values! referred to aero in the lowest rotational 
level of tb-? v ■ 0 state, are tabulated next. The Raaan shift is calculated from 
the difference in these values. The frequency shifts, and wavelength for the 
assumed pump frequency are shown next for the v ■ 0 level for rotational transi- 
tions, and v ■ 0 ■> 1 vibrational-rotational transitions. Lastly the populations 
and Raman intensities for the various Raaan branches are tabulated. The popula- 
tions are in units of cm“ 3 , and the intensities are ca~ 1 sR~ 1 . Energies and 
frequencies are in cm” * , and wavelengths are in Angstroas. 
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21724.63 

24615.56 


13684.37 

17082.43 

20204.04 

23048.17 

25600.70 


15197.73 

13490,51 

21494,29 

24204.59 

26616.93 


16792.43 

19967.04 

22836.49 

25395.19 

27637.57 


18481.26 

21522.14 

24238.55 

26623.74 

20670.97 


20291,67 

23170.01 

25711.07 

27896.86 

29719.43 


22216.37 

24929.50 

27263.63 

27224.49 

30707.84 


24314.00 

26324.66 

28930.29 

30620.31 

31CC4.09 


26609.86 

23805.40 

30720.19 

32102.21 

33017.45 


29146.61 

31148.18 

3266C.00 

33692.84 

34208.90 


31975.08 

33656.60 

34307.93 

35420.07 

35472.01 


35155.05 

36462.24 

37108.62 

37317.59 

36832.54 


38756.06 

39625.31 

39054.51 

39425.41 

30319,76 



R»»»n 

rraauancv Shift* 




J Rot-0 

Rot-S 

Vib-Q 

Vib-Q 

Vib-S 

0 

0.00 

354.37 

0.00 

4161*18 

4497.05 

1 

0.00 

507.02 

0.00 

4155.26 

4712.91 

*“» 

-354.37 

814.38 

3806.81 

4143.47 

4917.03 

3 

-587,02 

1034.57 

3563.24 

4125.89 

5100.46 

4 

-814.38 

1245.95 

3329.09 

4102.65 

5205*72 

5 

-1034.57 

1447.16 

3071.32 

4073.89 

5447.60 

# 

o 

-1245.95 

1637.24 

2856.71 

4039.77 

5593.25 

7 

-1447.16 

1815.65 

2626.73 

4000.44 

5722.19 

9 

-1637,24 

1982.33 

2402.53 

3956.01 

5834.36 

9 

-1815.65 

2137.83 

2134.80 

3706.55 

5930.15 

10 

-1982.33 

2283.31 

1973.68 

3852.03 

6010.43 

1! 

-2137.83 

2420.63 

1763.72 

3772.32 

6076.60 

12 

-2283.31 

2552.45 

1568.72 

3727.13 

6130.57 

13 

-2420.63 

2602.22 

1371.69 

3655.97 

6174.82 

14 

-2552.45 

2814.34 

1174.66 

3578.13 

6212.40 

15 

-2682.22 

2754.16 

973.74 

3492.59 

6246.94 




A-2 




M U A U M KO'OOMM^UMHOyjO^JO'WJ.WMM 


14 


•2814.34 


3108.05 


743.70 


3398.05 


4282.67 




17 

-2954 « 14 

3283.53 

538.44 

3292.78 

6324. 

41 

18 

-3108.05 

3489.24 

290.00 

3174*61 

6377. 

59 

19 

-3283.53 

3735.11 

9.25 

3040.88 

6448. 

24 

20 

-3489*24 

4032.33 

-314.63 

2880*34 

6542. 

99 

•’1 

-3735.11 

4373.49 

694.23 

2713.13 

6667. 

04 

22 

-4032.33 

4832.61 

• 1143.?? 

2510.65 

6834. 

18 

23 

-4393.49 

5365.23 

1680*36 

2275.55 

7046. 

74 

24 

-4832.41 

6008*44 

■2321.96 

2001.57 

7315. 

63 

25 

-5345.23 

4780.97 

-3039.68 

1601.52 

7650. 

22 



Ramon 

Uavalanathi 





Rotational 


Vibrational 


J 

Lamda-0 

Lamda-S 

Lamda-0 

Lamda-Q 

Laaida 

-8 

0 

0.00 

4765.83 

0.00 

6123.39 

6252. 

28 

1 

0.00 

5023.87 

0.00 

6121*17 

6337. 

50 

2 

4797.00 

5081.92 

5973.34 

6116.76 

6420* 

56 

3 

4744.05 

5139.42 

5908.05 

6110*19 

6500. 

46 

A 

4493.43 

5195.87 

5326.52 

6101*52 

6576. 

23 

w 

4445.42 

5250.77 

5746.90 

6090.84 

6646. 

99 

4 

4600.25 

5303.70 

5670.45 

6073*20 

6711. 

97 

? 

4558.06 

5354.36 

5597.45 

6063.71 

6770.57 

e 

4518.91 

5402.58 

5528.08 

6047.42 

6822. 

38 

? 

4482.77 

5448.35 

5462.33 

6029.38 

6867. 

26 

10 

4449.52 

5491.88 

5400*06 

6009.63 

6905. 

33 

i X 

4418.95 

5533.62 

5340.94 

5980.14 

6937. 

03 

12 

4390.72 

5574.27 

5204.50 

5964.85 

6963. 

10 

13 

4364.40 

5614.89 

5230.04 

5939.64 

6984. 

62 

1 4 

4339.44 

5656.86 

5176.70 

5912.31 

7003. 

00 

13 

4315.14 

5701.95 

5123.41 

5882.56 

7019. 

98 

14 

4290.68 

5752.43 

5063.88 

5350.02 

7037. 

63 

1 7 

4265.09 

5811.09 

5011*64 

5814.22 

7058.36 

; ?2; 

4237.28 

5881.40 

4750.00 

5774.54 

7084. 

96 

:? 

4206.00 

5967.69 

4882. 16 

5730.29 

7120. 

60 

20 

4169.93 

6075.46 

4306*16 

5630.64 

7163. 

97 

21 

4127.61 

6211.75 

4720.05 

5624,65 

7234. 

34 

• *» 

4077.58 

6305.94 

4621.93 

5561.32 

7321. 

61 

'3 

4018.41 

6610.79 

4510.12 

5409.54 

7437. 

57 


3948.73 

6704.38 

4303.28 

5400.20 

753?. 

34 

-> «r 

3867.3? 

7293.40 

4240.50 

5316.10 

7787. 

09 


Praaaura 

* 5500.00000 raia 




Tern* 

* 300 Qrot 

> 1.8882 

Qvib ■ 1 

.0000 





Rotational 

Vibrational 


J 

N < J ) 

I (0-branch) 

Z <S branch) 

Z(0 branch) Z(Q -branch) 

Z (8-branch 

0 

1.212E+21 

0 ♦ 000E400 

4.023E-0? 

0 . 000E4O0 

1.343E 0? 

1.333E-10 

1 

6. 086E+21 

0. 000E400 

1 . 409E-0B 

0.000E400 

8. 073E-09 

3.462C -10 

*> 

1 • 052E421 

9.613E- 10 

2 . 042E-09 

3.506E-U 

1.377E-0? 

4.397E-11 

3 

7.670E+20 

9.572E-10 

1 . 3S9E-09 

3.737E-11 

9.992C-10 

2.532E-11 

4 

3.183E419 

4.717E- 11 

5.336E -11 

1.949E-12 

4.129E 11 

8.497E-13 

J 

6. 304E+13 

1 . 06SE- 1 1 

1 . 031E-1 1 

4.606E-13 

8* 137E-12 

1.372E-13 

TiBil* 

■ 400 Qrot 

* 2.4503 

Qvib * 1 

.0000 





Rotational 

Vibrational 


J 

N< J) 

I (0-branch> 

Z (G -branch) 

Z(O-branch) : 

1(0 -branch) 

Z(6 -branch 

0 

6.983E+20 

0. 000E400 

2.778C-09 

0. 000E400 

C.890E-10 

7.6C1E- 11 

1 

4.056E+21 

0.000E400 

9.392E -09 

0.000E400 

5.381E-09 

2.306E-10 

A 

9.389E+20 

8.S79E-10 

1 • 822E-09 

3.129C- 11 

1 . 229E-09 

3.926E-11 

3 

1.060E421 

1 . 323E* 0 / 

1 . 377E-09 

5. 166E-11 

1.382E-0? 

3.501E-11 

4 

7.889E419 

1 » 169E-10 

1 . 327E-10 

4.C30E-12 

1 • 023E-10 

2.106E-12 

5 

3.241E419 

5.476E* 11 

5.301E-11 

2.360E-12 

4 • 163E-1 1 

7.053E-13 

6 

9.23- ?417 

1.750E-12 

1 • 493E-12 7.C42E-14 

A- 3 

1 . 1D5E-12 

1 • 638E-14 


'4J1? 


OfflOMM. PAGE It 
OP POOR QUALITY 

Taar a 300 Qrot • 3*0200 Qvib • 1*0000 

Rotational Vibrational 


J 

N< J) 

KO-branch) 

X(8-branch) 

Z(O-branch) 

KQ-branch) 

K8-braneh 

0 

4.535E420 

0aQ00E400 

1.804E-09 

0 *0008400 

3.7748*10 

4.9698-11 

1 

2.876E421 

O.OOOE400 

6.650E 09 

O.OOOEfOO 

3*0158-09 

1.636E-10 

n 

▲ 

7 • 9308420 

7.246E-10 

1.5398*09 

2*6438-11 

1.038E- 09 

3*3168-11 

3 

1.165E421 

1 .4538-09 

2.0648-0? 

5.674E-U 

1*5178 -09 

3.8458-11 

4 

1.230E420 

1.823E-10 

2.069E-10 

7*3308-12 

1.595E-10 

3*2038-12 

3 

7 • 8288419 

1.3238 10 

1 *280E -10 

5.719E-12 

1*0108-10 

1*7038-12 

6 

3.771E418 

7# 149E-12 

6.098E-12 

3.203E-13 

4.039E-12 

6.690E-14 

7 

1 • 125E418 

2.3768- 12 

1.Q22E -12 

1.092E-13 

1.433E12 

1*6198-14 

T aar 

■ 600 Qrot 

• 3*5981 

Qvib ■ 

1*0000 





Rotational 

Vibrational 


J 

N( J) 

I(O-branch) 

I (8-branch) 

I(O-branch) 

KQ-branch) 

Z(8* branch 

0 

3.180E+20 

0.0008100 

1.265E-09 

0*0008400 

4.049E-10 

3.499E-11 

1 

2 » 138E421 

Oa 000E400 

4.9S0E-09 

0.000E400 

2.8368*09 

1*2168-10 

2 

6.625E420 

6.054E* 10 

1.286E 0? 

2*2008-11 

0.671E 10 

2.770E-11 

3 

1 . 1S9E421 

1.447E-09 

2.055E-09 

5.6488-11 

1.510E-09 

3.0278 *11 

A 

1.5468420 

2.2918-10 

2.601E-10 

9.466E-12 

2.0068-10 

4.1278-12 

*# 

1 a 31BE420 

2.226E-10 

2. 155E-10 

9.6278-12 

1.7018 10 

2*0678-12 

6 

9 • 010E418 

1 . 700E-11 

1.457E-11 

7*6518-13 

1.1568*11 

1.5908-13 

7 

4 a 043E418 

8.540E-12 

6 aSSlE-12 

3*9258 13 

5*1518-12 

5.618E* 14 

8 

1 a 479E417 

3.470E* 13 

2.420E-13 

1.621E-14 

1.8698-13 

1.7108-15 

TVBiA- 

■ 700 Qrot 

■ 4*1683 

Qvib * 

1.0002 





Rotational 

Vibrational 


J 

N( J) 

I (0-branch) 

KS-branch) 

I(O-branch) 

1(0 branch) 

I (S-r.-anch 

0 

2.353E420 

0.000E 1-00 

9 » 361E-10 

0*0008400 

2.9968 10 

2*5888-11 

1 

1 a 649E421 

O.OOOE400 

3 .8108-09 

0.0008400 

2* 1C7E-09 

9.3C1E 11 

2 

5.SS4E420 

5.075E 10 

1.07SE-09 

1 *0518-11 

7.2698-10 

2.3238-11 

J 

1 a 101E4 21 

1 * 374E-09 

1.932E-09 

5.3668-11 

1.435E-09 

3.636E-U 

4 

1.735E420 

2.572E-10 

2 .720E -10 

1 .0638 *11 

2.2518-10 

4.6338-12 

3 

1 a 822E420 

3.078E-10 

2.9798-10 

1*3318- 11 

2.3518-10 

3.9648-12 

6 

1 a 600841? 

3.0338 11 

2.587E-11 

1.359E-12 

2.053E-U 

2.8338-13 

•a 

4 

9.612E418 

2 a 030E-1 1 

1.557E-U 

9*3318-13 

1.2258-11 

1 *3838-13 

s 

4.908E417 

1 « 152E-12 

3.057E 13 

5.331E-14 

6.201E-13 

5.6748-15 

T fe Ha ► 

• 800 Qrot 

• 4.7389 

Qvib - 

1.0006 





Rotational 

Vibrational 


J 

N( J) 

I <0-branch> 

I (8-branch> 

KO-branch) 

Z (Q-branch) 

Z (8-branch 

0 

1 a 810E420 

0*000Ei00 

7.202E 10 

0.0008400 

2.3058-10 

1.991E-11 

1 

1 a 309E421 

O.0OOE4O0 

3.030E-09 

0*0008400 

1.7368-09 

7.447E-11 

2 

4 a 693E420 

4.289E 10 

9. 109E-10 

1.5*48-11 

6.143E-10 

1 . 963E-1 1 

3 

1.022E421 

1 » 276E-09 

1.812C-09 

4* V60E-11 

1*3328-09 

3.3748-1 1 

4 

1 a 825E420 

2.705E 10 

3.071E* 10 

1.1188-11 

2.360E -10 

4.872E-12 

5 

2.240E420 

3.785E-10 

3.664E-10 

1.6378-11 

2*8928-10 

4.0758 12 

4 

2a 373E419 

4.499E* 11 

3.830E-11 

2.016E -12 

3.045E 11 

4.210E-13 


1 a 775E419 

3.749E-11 

2.876E-11 

1 • 723E-12 

2.2618 11 

2.5548-13 

8 

1 • 1648416 

2.7318-12 

1.911E -12 

1*2768-13 

1.4718-12 

1.3468-14 

? 

5.443E417 

1.417E-12 

9 a 138E-13 

6*6738*14 

6.810E-13 

4.9028-15 

Temr- 

« 900 Qrot ■ 5a3096 Qvib ■ 

1.0013 





Rotational 

Vibrational 


J 

N( J) 

KO-branch) 

Z (3-branch) 

KO-branch) 

KQ-branch) 

KS-branch 

0 

1 a 435E420 

0.0008(00 

5.709E-10 

0*0008400 

1.Q27E-10 

1.5798- 11 

1 

1.063E421 

0.000E )00 

2.462E-09 

0*0008400 

1*4108-09 

6.0498-11 

2 

4 . 002E+20 

3.657E-10 

7.76CE-10 

1.3348-11 

5.2388-10 

1.6748-11 

3 

9 a 377E420 

1 a 170E-09 

1.662E-09 

4.5638-11 

1.222E-09 

3.095E-11 

4 

1 a 845E420 

2.735E-10 

3* 105E-10 

1 a 1308-11 

2.3948-10 

4.9268-12 

3 

2 * 558E420 

4.3228- 10 

4*1 84E 10 

1 .3698-11 

3.3028-10 

5.5668-12 

6 

3.136E419 

5* 944E-1 1 

5.070E-11 . , 2.663E-12 
A-4 

4.0238-11 

5.563E-13 


7 

2.780E419 

S.873E 11 

4.S05E-11 

2.699E-12 

3.542E-11 

4.001E-13 

8 

2.215E418 

5* 197E-12 

3.636E-12 

2.420E-13 

2.7V8E-12 

2.560E 14 

9 

1 * 269E4 18 

3.356E- 12 

2.164E-12 

1 . 501E-13 

1 .613E-12 

1.180E-14 

10 

6.786E416 

1 « 960E-13 

1.175E-13 

9.232E-13 

0 . 396E-14 

4 * 819E-16 

T«.*r 

■ 1000 Qrot 

• 5.8802 

Qvlb ■ 

1.0025 





Rotational 

Vibrational 


J 

N( J) 

I <0-branch> 

KG-branch) 

Z(O-branch) 

Z(Q-branch) 

X (8-branch 

0 

1 * 165E420 

O.O0OE4O0 

4 • 634E-10 

0. OOOEiOO 

1.403E-10 

1.281E-11 

1 

8 » 799E420 

0 . 000E400 

2 . 037E-09 

O.0OOE4OO 

1 • 167E-07 

5.006E-11 

* 

3 ♦ 444E+20 

3. 147E-10 

4.484C-10 

1 . 148E-11 

4.507E-10 

1.440E-11 

3 

6* 554E420 

1.047E-07 

1.516E-09 

4.167E -11 

1 . 1 14E-09 

2.324E-11 

4 

1 . 820E420 

2.697E-10 

3 • 062E-10 

1 . 114C-11 

2 .360E-10 

4.857E-12 

5 

2« 780E420 

4.677E-S0 

4.546E-10 

2 .031E-11 

3.50QE-1O 

6.047E-12 

6 

3.830E419 

7.260E-11 

6. 192C-11 

3.252E-12 

4.914E-11 

6. 7742-13 

7 

3.891E419 

8.217E- 11 

6 . 304E-1 1 

3 .777E-12 

4.757E 11 

5.597E-13 

3 

3 . 622E418 

8 « 498E-12 

5.945E-12 

3.970E-13 

4.576E-12 

4.187E 14 

9 

2.511E413 

6.530E-12 

4. 216E-12 

3.077E-13 

3.142E12 

2.297E-14 

10 

1 • 606E417 

4 ♦ 639E-13 

2.7C0E-13 

2. 183E-14 

1.907E-1? 

1.141E- IS 


* 1100 Qrot 

« 6.4511 

Qvlb » 

1.0044 





Rotational 

Vibrational 


J 

N< J) 

I <0 -branch) 

I <S-branch> 

X (0 branch) 

X (Q-branch) 

I (8 -branch 

C 

9.634E419 

0.000E400 

3.033E- 10 

0. OOOEiOO 

1.227E- 10 

1 * 060E-1 1 

1 

7.395E420 

0 • OOOEiOO 

1.712E-0? 

0. OOOEiOO 

7.009E- 10 

4.207E-11 

n 

*» 

2 ♦ 788E420 

2.730E-10 

5.779E-10 

7.758E-12 

3.91 IE-10 

1.247E 11 

3 

7.784E420 

7.714E- 10 

1.300E 07 

3.792E-11 

1.014E-07 

2.567E-11 

4 

1 • 76SE420 

2 . 613E-10 

2. 9692-10 

1 . 080E-11 

2.209E- 10 

4 . 71 IE 12 

t# 

2 ♦ 919E420 

4.9332- 10 

4.775E-10 

2.133E -11 

3.763E-10 

6.353E-12 

W 

4 « 425E419 

8. 388E-11 

7.155E-11 

3.75BE -12 

5.67CE-11 

7.050E- 13 

7 

5.026E+19 

1.062E 10 

0.143E It 

4.079E-12 

6.403E 11 

7.232E-13 

3 

3.313E418 

1.247E-11 

0 . 722 r -12 

5 • 625E-13 

6.713E- 12 

6 • 1432-14 

9 

4.252E410 

1.107E 11 

7.137E12 

5 .2132-13 

5.3202-12 

3 .C92E- 1 4 

10 

3. 138E417 

9. 209E-13 

5.519E-13 

4.338E 14 

3.945E-13 

2.2642 15 

11 

1 . 816E+17 

5.82GE- 13 

3.267E 13 

2.726E-14 

2 . 221E-13 

7.737E-16 

Tam* 

■ 1200 Qrot 

• 7.0219 

Qvib » 

1.0069 





Rotatior al 

Vibrational 


J 

N( J) 

I (0-br#nch> 

L (S branch) 

I (0-branch) 

X (Q branch) 

I (S branch 

0 

8. 093E4 19 

0. OOOEiOO 

: . 220E -10 

0. OOOEiOO 

1.030E-10 

3.903E-12 

1 

& 4 29SE420 

0. OOOEfOO 

! .457E-09 

0. OOOEiOO 

0 . 3S0E-10 

3.5812 11 

2 

2.612E420 

2.387E- 10 

5.067E-10 

0.705E-12 

3.413E-10 

1 . 092E-11 


7 • 0S1E+20 

8.836E-10 

1 .255E-09 

3.449E-11 

7.225E- 10 

2.337E-11 

4 

1.693E420 

2.507E- 10 

2.040E-10 

1.036E-11 

2. 196E-10 

4.517E-12 

IT 

"J 

2.992E420 

5. 056E-10 

4.874E-10 

2 . 106E-11 

3.862E-10 

6. SHE 12 

6 

4.9UE419 

9.310E- 11 

7 • 941E-1 1 

4 » 171E-12 

6.302E-11 

3.713E-13 

7 

6 . 121E419 

1 » 293E-10 

9.918E-11 

5.742E- 12 

7.779E 11 

C.808E-13 

a 

7 . 196E418 

1.633E- 11 

1 . 131E1 1 

7.003E-13 

9.091E -12 

3.319E-14 

9 

6 . 488E418 

1.689E-11 

1 . 009E-11 

7.955E-13 

8.U8E-12 

5.940E- 14 

to 

5.556E417 

1.605E-12 

9.617E-13 

7.553E-14 

6.073E-13 

3.946C-15 

11 

3 . 667E417 

1 . 175E-12 

6 • 587E-13 

5 . 497E-14 

4.478E-13 

1.763E 15 

T®»r 

■ 1300 Qrot 

» 7.5728 

Qvib ■ 

1*0101 





Rotational 

Vibrational 


J 

N( J) 

I (0-branch) 

1(8 -branch) 

X (0-branch) 

X (Q-branch) 

I (3 -branch 

0 

6.836E+1? 

0. OOOEfOO 

2.740E-10 

0. OOOEIOO 

0.768E-11 

7. 5762-12 

1 

5 • 4 17E420 

0. OOOEiOO 

1 • 254E-09 

0. OOOEiOO 

7.136E 10 

3.082E-11 

*s 

A. 

2 • 299E420 

2 • 100E-10 

4*. jIE-10 

7.661E-12 

3.008E-10 

7.6122-12 

3 

6 . 445E420 

3.043E- 10 

1.142E-09 

3.140E-U 

3.396E 10 

2 . 127E-11 

« 

•4 

1.612E420 

2.300C- 10 

2 . 712E-10 

9.867E- 12 

2.090E-10 

4.302E -12 

5 

3 « 013E420 

5.071E- 10 

4.720E-10 

2 .2022-1 1 

3.867E -10 

4.557E-12 

iJ 

5* 290E419 

1 » 003E-10 

G.554E-11 

4 • 493E-12 

6 . 7C8E- 1 1 

9.30SE- 13 

9 

/ 

7 * 132E419 

1.506E- 10 

1.156E-10 4.724E-12 

A-5 

7.0B7E-1 1 

1 . 026E-12 


ORIGINAL PAGE IS 
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3 

?. 172E+18 

2* 152E-1 1 

1 « 506E-1 1 

1 . 005C-12 

1. 159E-11 

1 . 060013 

9 

9 • 150E+18 

2.332E- 11 

1.336E-11 

1 . 122E- 12 

1.14SE 11 

0.376E-14 

10 

8 * 763E+17 

2.332E-12 

1 .517012 

1 . 192013 

1.004C-12 

6.225E-15 

tl 

6. 545E+17 

2.093E- 12 

1.176E-12 

9.312E-14 

7.994E- 13 

3.504E-1S 

« 

• A* 

4.710E+16 

1.676E-13 

8.829E-14 

7.737E 15 

5.670E-14 

I.C27E -16 

Ten* 

* 140C Qrot 

* 8.1635 

Qvib ■ 

1.0141 





Rotational 

Vibrational 


J 

N( J) 

I(0-brar»ch) 

I<S -branch) 

KO-branch) 

X (Q branch) 

I (S-branch 

0 

5.924E+19 

0.000E+00 

2.357E-10 

0.0000 00 

7.543E-11 

6.517E-12 

1 

4.705E+20 

0.000E400 

1 .089E-09 

0.000E+00 

6.241E-10 

2.677E-11 

n 

*• 

2. 035E+20 

1 • 860E-10 

3.730E-10 

6. 783E-12 

2.6640-10 

8.510E-12 

3 

5. 874E+20 

7.330C- 10 

1 * 041E-09 

2.361E-11 

7.652E -10 

1.939E-11 

4 

1 ♦ S26E+20 

2.262E-10 

2.568E-10 

9.345E-12 

1.980E-10 

4.074C 12 

5 

2.994E+20 

5.060E 10 

4 * 890E- 10 

2. 180E-11 

3.365E-10 

6.516E- 12 


5.370E+19 

1 .0S6E-10 

9.006E-11 

4.730E-12 

7.146E- 11 

9.881E-13 

7 

8.032E+1? 

1.697E-10 

1.301E-10 

7.777E-12 

1.023E-10 

1 . 156E-12 

3 

1 • 116E+19 

2. 618E-11 

1 . 831E-1 1 

1 .223012 

1. 409E-1 1 

1 .270E-13 

9 

1.214E+19 

3.1S9E-11 

2.037E -11 

1 • 433E-12 

1 .S13E-11 

1.11 IE-13 

10 

1.280E+18 

3 . 697E-12 

2 . 216E-12 

1.741E 13 

1 .5030-12 

9.070E-15 

11 

1 ♦ 362E+1B 

3.405E- 12 

1.907E-12 

1.593E-13 

1.29CE -12 

5.6S9E-15 

1 -4 

3 • 53 IE 4 16 

3.034E-13 

1 . 609E-13 

1.410E-14 

1.033E-13 

3.32SE 16 

Tan* 

* 1500 Qrot 

« 3.7347 

Qvib * 

1.0189 





Rotational 

Vibrational 


J 

N( J) 

I (0-branch) 

I (S-branch) 

I ( O-branch ) 

I (Q -branch) 

I (5- branch 

0 

5 • 143E+19 

O.OOOE+OO 

2.046E- 10 

O.OOOE+OO 

’..549E-11 

5.65CC 12 

1 

4 • 119E+20 

O.OOOEIOO 

9.537E-10 

O.OOOE+OO 

5.464E-10 

2.344C 11 

*> 

1 . 812E+20 

1 • 656E-10 

3. 516E-10 

6. 038E-12 

2.371E-10 

7 . 3760-12 

3 

5.361E+20 

6.670E- 10 

7.302E-10 

2.612E 11 

6.985E-10 

1.770E-11 

> 
■ / 

1 » 440E+20 

2 « 135E-10 

2.424E-10 

0.C19E-12 

1.869E-10 

3.S45E-12 

u 

2.946E+20 

4.973E 10 

4.C19E-10 

2.153E-11 

3.C03E- 10 

6.411E-12 


5.761E+19 

1.092E-10 

9 • 316E-11 

4,0930 12 

7.372C- 11 

1 .0220 12 

-? 

4 

B.8C8E+19 

1.860E-10 

1.427E 10 

C.550E-12 

1 . 122E-10 

1.2670-13 

.J 

1 * 303E+19 

3.063E-11 

2 . 147C-1 1 

1 . 4340-12 

1.6520 11 

1.512E-13 

9 

1 ♦ 533E+ 1 7 

3.9720 11 

2.574E 11 

1 .030012 

1.91CE- 11 

1.4040-13 

10 

1.759E+18 

5 • 077E-12 

3 » 043E-12 

2 • 27 IE 13 

2.175E-12 

1. 2430 14 

• 4 

A A. 

1 .599E+18 

3. 126E 12 

2.C73E- 12 

2 • 373E- 13 

1.953E 12 

C.563E-15 

■•» 

1 * 428E+17 

5.081E-13 

2 .676013 

2.345E-14 

1.7190 13 

5.5370 16 

•■ T 

A o 

1 .014E+17 

4.0112 13 

1.771E 13 

I.C17E-14 

1.201E-13 

2.63CE- 16 

T -r 

= 1600 Qrot 

= 7.3055 

Qvib * 

1.0244 





Rotational 

Vibrational 


i 

•J 

NCJ) 

I (9-branch> 

I (S-branch) 

I (P- branch) 

I (Q-branch> 

I (S- branch 

0 

4.502E+19 

0.000000 

1.791E-10 

O.OOOE+OO 

5.732E-11 

4.952E -12 


3.631E+20 

O.OOOE+OO 

3. 40DE-10 

O.OOOE+OO 

4.817C 10 

2.066E-11 

4. 

1 .621E+2G 

1.431E- 10 

3.146E -10 

S.402E-12 

2.121E 10 

6.778E-12 

3 

4 . 902E+20 

6.118E-10 

8 . 689C-10 

2. 388E-11 

6.387E-10 

1.618C-11 

4 

1 • 356E+20 

2.010E- 10 

2.262E -10 

0.303E-12 

1 . 759E-10 

3.620E-12 

’J 

2 .877E+20 

4.861E-10 

4 . 705E-10 

2.102011 

3.713C-10 

6. 260E-12 

6 

5.C77E+19 

1 . 114E-10 

7.503E 11 

4. 971E-12 

7.541E-11 

1 .0430-12 

f 

/ 

9.456E+19 

1.997E-10 

1 • 532E-10 

9 . 100012 

1 * 205E-10 

1.3610 12 

B 

1 .48BE+19 

3.492E- 11 

2.443E -tl 

1.632E-12 

1.330E- 11 

1.721E-13 


1 ♦ 864E + 19 

4 . 832E- 1 1 

3.129E-11 

2.285E-12 

2.332E-11 

1 • 706E-13 

i ■; 

2 .298E+18 

6.633E- 12 

3.97CE-12 

3 . 127E-13 

2.C43E-12 

1.632E-14 


2.246E+13 

7.262E-12 

4 ♦ 070E-12 

3.377C-13 

2.767E-12 

1.2130-14 

« o 

2 . 208E+ 17 

7.857E 13 

4. 13CE -13 

3.627E-14 

2.653E-13 

6. 562E-16 

1 5 

1 . 724E + 17 

6.819E-13 

3.385E-13 

3.0C9C- 14 

2. 042C-13 

4.370E-16 


»»*«»• 
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